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ABSTRACT 

The conversion of the globular cluster luminosity function (GCLF, dN/d\ogL) to the globular 
cluster mass function (GCMF, dN/d\ogM) is addressed. Dissolving globular clusters (GCs) become 
preferentially depleted in low-mass stars, which have a high mass-to-light ratio. This has been shown 
to result in a mass-to-light ratio (M/L) that increases with GC luminosity or mass, because more 
massive GCs have lost a smaller fraction of their stars than low-mass GCs. Using GC models, we 
study the influence of the luminosity dependency of M/L on the inferred GCMF. The observed GCLF 
is consistent with a powerlaw or Schechter type GC initial mass function in combination with a cluster 
mass-dependent mass loss rate. Below the peak, the logarithmic slope of the GCMF is shallower than 
that of the GCLF (0.7 versus 1.0), whereas the peak mass is 0.1 — 0.3 dex lower when accounting for 
the variability of M/L than in the case where a constant AI/L is adopted. 

Subject headings: Galaxy: kinematics and dynamics — (Galaxy:) globular clusters: general — 
galaxies; kinematics and dynamics — galaxies: star clusters — stellar dynamics 



1. INTRODUCTION 

The present-day globular cluster mass function 
(GCMF, dN/dlogM) is derived from the globular clus- 
ter luminosity function (GCLF, dN/dlogL) by as- 
suming a constant mass-to-light ratio for all globular 
clusters (e.g [Fall fc Zhang 200 it iVesperini et al.l [200l 
[Jordan et al.] 120071 : EcLaughlin i^Fall 12008!) . The re- 
sulting GCMF is strongly depicted in low-mass globular 
clusters (GCs) with respect to the mass distribution of 
young star clusters, which is well-described by a pow- 
erlaw with index —2 in various environments down to 
a few 100 Mq. This has led to a number of pioneer- 
ing studies explaining its shape by cluster evaporation 
at a cluster mass-independent mass loss rate (equiva- 
lent t o a disruption time ^dis oc M, e.g., iFall fc Zhani 
120011: lVesperinill2001h acting on a powerlaw or lSchechtei 
' 1976f l cluster initial mass function (CIMF, e.g. , 



Harris fc Pudritd 119941: [Mc Laughlin fc Pudritz 





Elmegreen fc Efrcmovl 119971: IBurkert fc SmithI 
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Gieles et al.. 2006 ). 



Although the observed peaked shape of the GCMF 
is reproduced in the above studies, the underlying 
assumptions arc not entirely satisfactory because a 
cluster mass-dependent mass loss rate (equivalent to 
^dis oc M'^ with 7 ~ 0.7, se e Eq. [J) is found in 
theory (e. g., Q jaumgardt 2001; iBaumgardt fc Makinol 



iory 

2003HGieles fc Baumgardt,2008i) and observations (e.g. , 
Lamers et al .1 120051: IGieles fc Bastianl[2008l : ILarsenl[2009l : 



Gielejl2009D . This arises from the non-linear scaling of 



the disruption time with the half-mass relaxation time 
(tdis oc ^)j which is caused by the non-zero escape time 
of stars with velocities above the escape velocity from a 
tidally hmited cluster (jFukushige fc He ggid [20001) . The 
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physical effect of a lower 7 is that the dissolution rate 
of low-mass clusters is slowed down relative to higher 
cluster masses and higher 7. 

The low-mass slope of a dissolution-dominated mass 
functio n like the GCMF i s always equal to the expo- 
nent 7 ([Fall fc Zhang|l200 II : [Lamers et al.ll2005( ). A mass- 
dependent mass loss rate conflicts with the observations, 
as it yields a higher number of low-mass GCs compared 
to cluster mass- independent mass loss (7 = 1). The dis- 
agreement between cluster mass-dependent mass loss and 
the observed sparse population of low-mass GCs is illus- 
trated in Fig. dja) . The slope of the modeled low-mass 
GCMF is ^ 0.7 for mass-dependent mass loss (7 = 0.7), 
whereas for cluster mass- independent mass loss (7 = 1) 
the slope is 1.0, in agreement with the observations. 
The peak (or 'turnover') masses also differ by ~ 0.3 dex. 
These differences show that a lower mass loss rate for 
low-mass GCs (7 = 0.7) yields a higher number of these 
relative to massive GCs than in the case of a constant 
mass loss rate (7 = 1). 

Recent studies show that the mass-to-light (M/L) ra- 
tios of GCs are not constant with lu minosity or mass 
(jReikuba et all 120071 : iKruiissenI [2008f) . contrary to the 
assumption of a constant M/L ratio ir i previous studies. 
This a grees with an earlier analysis bv lMandushev et "all 
(|1991h . who determined dynamical masses of Galactic 
GCs. These studies show that M/L increases with mass 
and luminosity because low-mass GCs are more strongly 
depleted in low-mass stars. This variation of M/L will 
affect the conversion of the GCLF to a mass function. 
Specifically, the smaller M/L ratios of low-mass GCs 
imply that the masses of low-mass clusters are overes- 
timated and consequently, that the low-mass end of the 
GCMF would be shallower than presently expected. The 
variability of M/L could therefore strongly affect the in- 
terpretation of the GCLF. 

We show that the relation between the GCLF and 
the GCMF is affected by low-mass star depletion, which 
arises from two-body relaxation (e.g., Meylan fc Heggii 
Il997f ). In Sect. [2] we discuss the influence of a luminosity- 
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Fig. 1.— The inferred GCMF of Galactic GCs I IHarrisI \TMEl . 
histogram) . Panel (a): GCMF derived using M/Ly = 3 (as in 
I Fall &: Zhand I2OOII ') . Overplotted is our model MF with a mass- 
dependent mass loss rate (solid line, see Sect.[2]l ado pting a disso- 
lution timescale to = 1.3 Myr. Using t^is = toM"' l ILamers et al.l 
120051") . for a 10^ Mq GC and 7 = 0.7 this corresponds to a dis- 
ruption time of tjis = 21 Gyr. The dashed line sh ows the model 
for a cluster mass-independent mass loss rate (as in ll'all Ik Zhand 
[2001D . Panel (b): GMCF derived from the GCLF using the lumi- 
nosity dependent M/Ly (see Fig. |4jl. The solid curve is the same 
as above while the dotted curves represent models for (from bottom 
to top) log (to /Myr) = log 1.3 + {-0.5, -0.25, 0.25}, corresponding 
to t(jis = 7 — 37 Gyr. Error bars are Itr Poissonian. 

dependent mass-to- light ratio on the inferred GCMF, 
and we model the GCLF in Sect. [3l alleviating the ob- 
servationally expensive need for accurate M/L ratios to 
derive the GCMF to allow for a comparison with theory. 
By including a mass-dependent mass loss rate and a vari- 
abl e M/L ratio, our mo del provides an improvement to 
the lFall fc Zha^ (pOOll ) model. 

2. IMPLICATIONS OF A LUMINOSITY-DEPENDENT M/L 

We model the evolution of star clusters in order to 
quantify the influence of the luminosity dependence of 
M/L on the relation betw een the GCLF and the GC MF. 
Our model, called SPACE (jKruiissen fc Lamer j|2008f ). in- 
cludes mass loss by stellar evolution and by evaporation. 
The mass loss by evaporation is pa rametrised with the 
simple relation (jLamers ei~all[20()5l) : 

/dM\ M _ M^-^ 



dissolution timescale which depends on the environment. 
We illustrate the effect of a variable M/L ratio by adopt- 
ing a unique value for to , which we assume to be the same 
for all clusters (a realistic spread in to is considered in 
Sect. [3]). We subsequently convert the observed LF of the 
sample of GCs to a MF by adopting the corresponding 
relation between Ly and M/Ly that is computed with 
SPACE (see Fig. g]). In Fig. [U ^b) we s how t he resulting 
MF for the 146 GCs from the iHarrisI (|1996f l catalogued 
Overplotted are the model MFs with diff erent values for 
to, adopting a metallicity Z = 0.0004, a iKroupal ()200lD 
stellar IMF and a Schechter CIMF with powerlaw index 
—2 and exponential truncation mass M* = 2.5 x 10^ Mq. 
As expected from Eq. [U the slope of the MF is indepen- 
dent of the dissolution timescale. 

By comparing panels (a) and (b) in Fig. [1] we see that 
the luminosity dependency of M/L gives rise to two ef- 
fects: (1) the slope at the low-mass end of the inferred 
GCMF drops to ^ 0.7, which is the expected value for 
mode ls with cluster a ma ss-dependent cluster mass loss 
rate (jLamers et al.l [20051 ) . and (2) the peak in the MF 
(the so-called turnover mass) shifts to a lower mass with 
^ 0.3 dex. About half this shift is due to the already 
high value of M/L = 3 adopted bv lFall fc Zhand (|2001). 
The slope of the GCMF at the low-mass end is different 
from the slope of the GCLF, and therefore also different 
from the GCMF slope (~ 1) that would be inferred from 
the GCLF when using a constant M/L ratio. 

3. MODELS OF THE GALACTIC GC SYSTEM 

In our above analysis we have assumed a single dis- 
solution timescale for the entire GC system. In reality 
there is a range of timescales on which the GCs dissolve. 
We now consider a more detailed Monte Carlo model of 
the Galactic GC system in which the dependency of the 
dynamical evolution of GCs on their orbits is included. 
Our aim is to directly model the GCLF, rather than to 
obtain it by converting the GCMF. 

The initial positions of the GCs with respect to the 
Milky Way are taken f rom the powerlaw-like density 
profile (see e.g., IML&Z hang 2001), Eq. 26) that arises 
from the isothermal sphe re, with an outward increase of 
the velocity anisotropy (jEddingtonI 119151 ). Our choice 
of parameters for the kinematic model are (1) an ini- 
tial anisotropy radius i?A = 1 kpc, (2) a circular ve- 
locity of the gravitational potential Vc = 220 km s~^, 
and (3) (V^/w)^ = 3.5, which determi nes the slope of 
the density profile (iFall fc Zhang! [200lb . with v denot- 
ing the radial velocity dispersion. The initial velocities 
of the GCs are assigned according to the corresponding 
velocity ellipsoid (Aguilar et al. 1988, Eq. 4), including a 
systemic rotation of V^ot — 60 km s~^. We do not claim 
that this is the correct kinematic model for the Milky 
Way, but we consider it an appropriate ansatz. The re- 
sulting dissolution timescales agree with the range that 
is expected from observations. For less anisotropy, the 
mean dissolution timescale of surviving clusters would 
be longer . The initial cluster masses are drawn from a 
iSchechteil (|1976D function with index —2 an d exponential 
truncation mass M^ = 3 x 10^ Mq (also see [Jordan et al.l 



Here 7 = 0.7 for clusters with a Ki ng parameter typi- 1 wc adopt the 2003 edition of the data, which is available online 

Cal to GCs of = 7 (,LamerS et al.L2009. ) and to is the at http://www.physics.mcmast6r.ca/~harris/inwgc.dat. 
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Fig. 2. — Histograms of the initial (dashed) and present-day 
(solid) distributions of dissolution timescales to, tot • 

120071: [Harris et all 120091 )^. We sampl e the metalli cities 
from their observed distribution in the lHarrii (|1996l ) cat- 
alogue. 

Th e GC orbits are i ntegra ted in the Galactic potential 
from iJohnston et al.l ()1995t ). consisting of a bulge, disc 
and halo. We adopt a 4*'^-order Runge-Kutta integra- 
tion scheme with a variable timestep, in which the an- 
gular momentum and energy are conserved within 10~^ 
during each timestep. To compute the evolution of a GC 
with a given initial mass and metallicity, we derive its 
instantaneous dissolution timescale from the orbital pa- 
rameters. Tidal evaporation due to two-body relaxation 
and disc shocks are the main d issolution mechanisms 
(iChernoff et al.lll986D . Following iBaumgardt k Makinol 
(|2003l ) for the dissolution timescale due to two-body re- 
laxation we write: 

where i?a is the apogalactic radius of the cluster orbit, 
Vc.a is the circular velocity of the gravitational poten- 
tial at i?a, and e is the orbital eccentricity. The dis- 
solution timescale due t o disc shocks is expressed as 
(jCnedin fc Ostrikeilll997l . Kruijssen et al., in prep.): 

io.disc = 7.35 Myr f-^) P2A-\x), (3) 

where 14,5 is the velocity in the ^-direction during disc 
crossing at 2: = in units of 10^ m s^^, P2 is the 
(radial) orbital period in units of 10^ Myr, (7m, -10 is 
the orbital maximum of the acceleration due to the 
disc — 9$disc/5z in units of 10"^" m s~^, and Av,{x) 
is t he IWeinberel (11 994a, b.cl adiabatic correction'^ (see 
also IGnedin fc Ostriker. .1997i) . For mathematical sim- 

^ This number is slightly larger than in Sect. [21 because the 
spread in dissolution timescales implies that surviving massive GCs 
on average have a smaller dissolution timescale than surviving low- 
mass GCs. We correct for the resulting deficiency of massive GCs 
by increasing M*. 

^ The parameter x im plicitly depends on the GC mass (e.g., 
IGnedin fc Ustrikerl liggTl) . We adopt 0.6 times the initial GC mass, 
in agreement with the average mass loss per Hubble time from 
[Kruijssen fc Mieska L20(M ). 
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Fig. 3. — Histograms of the observed (filled) and modeled (dot- 
ted) GCLFs. 

plicity, we assume a very weak mass-radius relation"* of 
rh oc M°-* (Larson 2004) in the derivation of Eq. [3l 

Over a timespan of 12 Gyr, the dissolution timescales 
due to tidal evaporation and disc shocks are computed for 
every orbital revolution, measured between subsequent 
passages of the apogalacticon. The dissolution timescale 
that describes the mass loss rate (see Eq. [ij due to both 
effects is determined by adding the averaged inverses of 
both timescales: 



^O.tot ^0,cvap ^O.disc 

The resulting initial and present-day distributions of 
io.tot are shown in Fig. [2l GCs with short dissolution 
timescales are easily destroyed, leading to the deple- 
tion of the quickly-dissolving end of the distribution (at 
low values of to). The surviving GCs have dissolution 
timescal es that are in excellent agreement with other 
studies (p<niiissen[[2008[ : [Kruiissen fc Mieskd [2009,1 . Al- 
though their mean galactocentric radius is a factor two 
smaller than that of the observed Galactic GC system, 
the slopes of both density profiles are comparable. 

The evolution of GC mass and photometry is computed 
with SPACE, using the setup discussed in Sect.^l In total 
507079 GCs are generated with initial masses M > 5 x 
10^ Mo, of which 2000 survive until t = 12 Gyr. The 
present-day mass and luminosity functions are scaled to 
match the observed numbers, of which the scale factor 
can be used to derive properties of the initial Galactic 
GC system (see below). The computed V-band GCLF 
is compared to the observed distribution in Fig. [31 The 
distributions are in satisfactory agreement, with a KS- 
test p-value of 0.02. At low luminosities there is a slight 
discrepancy, which could be caused by incompleteness 
due to obscuration by the Galactic bulge (Gieles et al., 
in prep.). 

The M/ Ly ratio s of our modeled GCs are comp a red to 
observations from [McLaughlin fc van der Mare][ ([2005l l 
in Fig. [31 If low-mass star depletion is neglected (panel 
(a)), the M/Lv ratios of the models are completely set by 
their metallicities and they agree poorly with the obser- 

* Compared to adopting a constant GC radius, this assumption 
effects a ~ 0.45 dex scatter of io.disc- Because log {to disc/to,tot) > 
0.75 for all surviving GCs, this does not affect our results. 
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Fig. 4. — Comparison of observed M / Ly ratios of Galactic GCs 
(thick dots, ^McLaughlin & van dcr Marel 2005) with our modeled 
M/Lv ratios (thin dots). Error bars are Icr. Dotted horizon- 
tal lines denote the constant M/Ly ratios that are expected if 
low-mass star depletion is neglected (from bottom to top Z = 
{0.0004,0.004,0.008,0.02}). Panel (a): thin dots represent mod- 
eled M/Ly ratios without low-mass star depletion. Panel (b): thin 
dots denote modeled M/Ly ratios including low-mass star deple- 
tion. The solid line represents the relation between M/Ly and Ly 
that was adopted in the simple model of Fig. mb) . 

vations. When including low-mass star depletion (panel 
(b)), the modeled M/Ly ratios are affected by dynamical 
evolution and are in good agreement with the observa- 
tions. The same approach can be used to explain the 
observations o f Cen A, M31 and the LMC compiled by 
iReikuba et al.l (|2007f l. which gives results that are con- 
sistent with_ou£ajial^[si^^ 

The lMcLaughlin k van der Marel (2005) sample is not 
representative of the entire Galactic GC population, as 
it lacks GCs that are much fainter than the turnover 
and represents central rather than global M/Ly ratios 
for certain GCs, only allow ing for a first-order com - 
parison (for a discussion, see IKruiissen fc Mie ske"2009') . 
The observed slopes of the low-mass ste l lar m ass func- 
tions of 20 GCs from iDe Marchi et al.l (|2007D provide 
an independent check. Their compilation exhibits a 
clear trend of mass function slope with GC luminosity. 
Splitting their sample at about the turnover luminosity 
(log(iv'/-^0) = 5.1), for a mass function n oc m^" the 
mean slopes in the stellar mass range m = 0.3 — 0.8 
are abright = 1-42 ± 0.10 and afaint = 0.56 ± 0.07 for the 
bright and faint GCs, respectively. Faint GCs are in- 
deed more depleted in low-mass stars than bright GCs, 
substantiating our model results. 
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Fig. 5. — Histograms of the mass distributions of the modeled 
GCs. Represented are the CIMF (upper solid) and the present-day 
GCMF (dotted, with Itr Poissonian error bars). The GCMF for a 
single dissolution timescale (Fig. [Ub)) is represented by the con- 
tinuous solid curve. The dashed line gives the GCMF that would 
be obtained from Fig.|3]if a constant M/ L ratio were adopted. The 
initial mass distribution of the surviving GCs is given by the lower 
solid line. 

The initial properties of the Galactic GC system are 
obtained by scaling the present-day number of modeled 
GCs to the observed number and applying the same scale 
factor to the CIMF. In Fig. [5] we show the CIMF, the 
modeled GCMF, the GCMF that would be obtained from 
Fig.[3]if a constant M/L ratio were adopted, and the ini- 
tial mass distribution of the surviving GCs. The modeled 
GCMF for a single dissolution timescale from Fig. [TJb) 
is overplotted for comparison, illustrating its acceptable 
agreement with our detailed model. The disagreement 
for GC masses < 10'^ M© is due to the use of logarith- 
mic timesteps in our models, causing some GCs to lose 
their last few 100 M© within a single timestep at large 
ages. For a lower mass limit of the CIMF of 5 x 10^ 
(10^) Mq, we find a surviving GC number fraction of 
3.9 (0.1) X 10"'^, with an initial total mass of about 1.1 
(1.8) X 10^ M© and a present-d ay mass of 2.8 x lO"^ Mq. 
If the stellar halo (~ 10^ Mq, iBell et al.|[2008[) is con- 
stituted by disrupted GCs and coeval star s (in spite of 
chemical analyses, e.g.. lGratton et aLll2000( ). our compa- 
rable initial total GC mass implies that either nearly all 
star formation occurred in clusters at the epoch of GC 
formation, or that most of these stars now constitute the 
Galactic bulge. 

4. DISCUSSION 

We have shown that the interpretation of the GCLF 
as a one-to-one representation of the GCMF is incor- 
rect. This follows from the M/L ratio decrease due to 
the low-mass star depletion that arises from two-body 
relaxation. There is no equivalence of the luminosity 
function and the mass function as both have intrinsi- 
cally different low-mass slopes (~ 1 and ~ 0.7, respec- 
tively). In addition, the turnover mass is overestimated 
by 0.1 — 0.3 dex if a one-to-one conversion from GCLF 
to GCMF is applied, depending on the adopted M/L 
ratio. We have shown that the present-day GCLF and 
GCMF arise from a cluster mass-dependent mass loss 
rate (^dis oc M" "^ and idis oc i^^)' starting from a 
Schechter-type CIMF. Therefore, neither is consistent 
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with a cluster ma ss-independent mass loss rate (e.g., 
iFall fc Zhani[200il ). The GCMF that is computed us- 
ing a spread in dissolution timescale to only marginally 
differs from that for a single, mean value of ^o- 

The low-mass slope of a dissolution-dominated mass 
function hke the GC MF is equal to 7 (see Eg. [B, inde- 
pend ent of the CIMF (|Fall fc Zhandl2001l : iLamers et all 
120051 ). For cluster mass-dependent mass loss (7 = 0.7) 
the GCMF that is inferred from the GCLF is accu- 
rately matched by the models (see Fig. [l{b)). To ver- 
ify whether this perhaps holds for all values of 7, we 
have al so considered cluste r mass-independent mass loss 
(7 = l. lFall fc Zhantdl2001[ ) and found that the luminos- 
ity dependency of M/L (see Fig. lU is steepened com- 
pared to cluster mass-dependent mass loss. The conver- 
sion of the observed GCLF to a GCMF then gives an 
inferred GCMF slope that is even lower 0.6), in bad 
agreement with the expected (~ 1) value. We conclude 
that the match between the models and the observations 
only exists for values of 7 « 0.7. Of course, the precise 
description of mass loss does not affect the fundamen- 
tal principle of low-mass star depletion due to two-body 
relaxation. The luminosity dependence of M/L flattens 
the inferred low-mass GCMF in any scenario. 

We have not yet considered the radial variation of 
the turnover luminosity LtOj which has been shown 
to be independent o f galactocentric radius in M87 
(jVesperini et al.l [20031 ). Our prescription for dynamical 
evolution in Sect. [3] yields a higher turnover luminosity 
near the galactic centre than at large distances. How- 
ever, our method is aimed at investigating the influence 
of a representative spread in dissolution timescales on 
our results, rather than making an exact model of the 
Galactic GC system. It should be emphasised that the 
difference between the GCLF and the GCMF persists, 
even though it remains to be explained why Lto ap- 
pears to be constant. It could be that the outer GCs 
dissolve more rapidly than expected. Potential expla- 



nations could be that G Cs on wide or bits originate from 
accreted dwarf galaxies (jPrieto fc Gne din 2008i) , or a dis- 
solu tion mechanism that has not yet been included (see 
also lKruijssen fc ACe skc 2009 i), like the dynamica l impli- 
cations of whit e dwarf ki cks (Fregcau et al. 2009"), stellar 
evolution (Vcs perini fc Zepf 2003; Vespe rini et al. 2009), 
or gas expulsion (|Baumgardt fc Kroupali2007l ). 

The results of this paper do not only apply to the Milky 
Way, but also to other galaxies. We see that the prop- 
erties of the inferred GCMF are affected by the mass 
and luminosity dependence of M/L that ensues from 
low-mass star depletion. It is advised for observational 
and theoretical studies to be cautious when comparing 
GCLFs and GCMFs. At present an observed GCMF can- 
not be accurately obtained, because for most observed 
GCs only photometric masses are determined (for which 
by definition a constant M/L ratio is assumed) instead 
of dynamical masses. Considering the intrinsically differ- 
ent shapes of the GCLF and GCMF, the presently most 
feasible way of comparing theory and observations would 
be if models of GC systems are aimed at explaining the 
GCLF rather than the mass distribution. 
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